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Abstract. Photometric observations in the R and / bands of the detached M-type double-lined eclipsing binary 
CU Cnc have been acquired and analysed. The photometric elements obtained from the analysis of the light 
curves have been combined with an existing spectroscopic solution to yield high-precision (errors S^2%) absolute 
dimensions: M A = 0.4333 ± 0.0017 Mq, M b = 0.3980 ± 0.0014 M 0 , R A = 0.4317 ± 0.0052 R 0 , and R B = 
0.3908 ± 0.0094 R 0 . The mean effective temperature of the system has been estimated to be T e g = 3140 ± 150 K 
by comparing multi-band photometry (optical and infrared) with synthetic colors computed from state-of-the-art 
model atmospheres. Additionally, we have been able to obtain an estimate for the age (~320 Myr) and chemical 
composition ([Fe/H] ~ 0.0) of the binary system through its membership of the Castor moving group. With all 
these observational constraints, we have carried out a critical test of recent stellar models for low-mass stars. 
The comparison reveals that most evolutionary models under estimate the radius of the stars by as much as 10%, 
thus confirming the trend observed by Torres & Ribas ( |2002| ) for YY Gem and V818 Tau. In the mass-absolute 
magnitude diagram, CU Cnc is observed to be dimmer than other stars of the same mass and this makes the 
comparison with stellar models not so compelling. After ruling out a number of different scenarios, the apparent 
faintness of CU Cnc can be explained if its components are some 10% cooler than similar-mass stars or if there 
is some source of circumstellar dust absorption. The latter could be a tantalizing indirect evidence for a coplanar 
(Vega-like) dusty disk around this relatively young M-type binary. 

Key words, binaries: eclipsing - binaries: spectroscopic - Stars: low-mass, brown dwarfs - Stars: fundamental 
parameters - Stars: individual (CU Cnc) - planetary systems: protoplanetary disks 


1. Introduction 


Low-mass stars constitute the most numerous stellar pop¬ 
ulation of the Galaxy. Yet, their physical properties are 
still poorly known. This is not without the efforts of 
many investigators that have been working intensively 
on the modelling of the atmosphere and interior struc¬ 
ture of sub-solar mass stars. The major problem that the 
held has to face is a conspicuous lack of suitable cali¬ 
brators. Unfortunately, very few main-sequence M stars 
have empirically-determined masses, radii, luminosities 
and temperatures. 

Eclipsing binaries have often served as valuable bench¬ 
marks for the validation of structure and evolution models 


(see, e.g., Andersen |l99l|) . Despite some unsolved prob- 
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lems at the high-mass end of the main sequence, the be¬ 
haviour of stars more massive than the Sun is well under¬ 
stood, partly thanks to the large number of eclipsing bi¬ 
naries with well-determined physical properties. The situ¬ 
ation is radically different at the low-mass end of the main 
sequence. Selection effects caused by the faintness of the 
stars and the often strong intrinsic variations due to mag¬ 
netic activity pose a serious challenge for the discovery and 
analysis of detached eclipsing binaries. Only three eclips¬ 
ing systems with M-type components have been identi¬ 
fied to date. The member of the Castor multiple system 
YY Gem has almost identical components of spectral type 
Ml Ve and mass ^0.6 Mq (Torres & Ribas 2002, here¬ 
after TR02). CM Dra holds the record for the lowest mass 
eclipsing binary known and has components of ~0.25 M 0 
and spectral type M4.5 Ve (Metcalfe et al. 199C). These 
were the only two M-type eclipsing binaries known for 
decades until Delfosse et al. ( 1999a , hereafter D99) re¬ 
cently reported the discovery of eclipses in CU Cnc. 
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CU Cnc (GJ 2069A, HIP 41824, a = 8^31 m 37!58, 
S = +19° 2 3'39'/5) is an 11.9-mag spectroscopic binary 
with M3.5 Ve components. Also, it has a fainter visual 
companion, CV Cnc or GJ 2069B, at an angular distance 
of ~12". GJ 2069B is itself a close binary star (Delfosse 
et al. 1999b| ) with a mean spectral type of M4 and ap¬ 
pears to form a physically bound quadruple system with 
CU Cnc. Interestingly, the masses found by D99 place CU 
Cnc nicely in between of YY Gem and CM Dra and make 
it a system of paramount importance to improve upon our 
knowledge of the physics of low-mass stars. 

Currently, the high-precision spectroscopic orbit of CU 
Cnc presented by D99 (errors in minimum mass below 
0.4%) is not matched by the photometric light curves, 
and the authors reported only sparse photoelectric obser¬ 
vations. The motivation of our study was to acquire and 
analyse new high-quality photometric observations with 
the aim of determining accurate masses and radii for the 
components of CU Cnc. With a careful estimation of the 
temperature, our goal is to compare the physical proper¬ 
ties of CU Cnc’s components with the predictions of the 
available low-mass stellar evolution models. 


2. Light curve observations 

The just recent discovery of the eclipsing nature of CU Cnc 
is likely to be a consequence of its faintness but also, very 
importantly, of its shallow and short eclipses. The photo¬ 
metric observations published by D99 indicate a system 
with an orbital period of 2.77 days, eclipses of only ~0.2 
mag in depth and about 2 hours in duration. However, the 
high scatter of these light curves and their poor phase cov¬ 
erage prevented an accurate determination of the physical 
properties of CU Cnc. 

We report here new photometry of CU Cnc in the 
Johnson R and I bands. Both the photometric accuracy 
(a few millimagnitudes) and the phase coverage (over 2100 
observations) are sufficient to guarantee a reliable deter¬ 
mination of the light curve parameters thus permitting 
a critical evaluation of stellar models. The observations 
were carried out with the Four College 0.8-m APT, which 
is equipped with a refrigerated Hamamutsu photoelectric 
detector and filters closely matching the standard Johnson 
system. Differential photometry was obtained in which 
HD 72093 (F8, V = 7.80) was employed as the com¬ 
parison star, while HD 72358 (F5, V = 8.31) served as 
the check star. Although it is common practice to employ 
comparison stars that match the spectral type of the vari¬ 
able star, we preferred to avoid G, K or M stars because 
of concerns with magnetic-related intrinsic variability. In 
contrast, late F stars are known to be photometrically 
very stable. Note that only precise R and I photometry 
could be obtained with our instrumental setup because of 
the faintness of CU Cnc at shorter wavelengths. The sys¬ 
tem is so red that at the I band, its brightness rises up to 
about I ~ 9 mag. No evidence of significant light varia¬ 
tions (down to the few milimagnitude level) was found for 
the comparison-check star sets. The observations were re¬ 


duced using photometric reduction programs at Villanova 
University (USA). Differential extinction corrections were 
applied, although these corrections were typically very 
small. 

There is one feature that makes our photometric 
dataset especially noteworthy: the ~2100 measurements 
were obtained over a period of only 110 days (from 
December 5, 1998 through March 25, 1999). Classified as 
a flare star, CU Cnc is expected to be magnetically active 
and display some level of brightness variability due to sur¬ 
face inhomogeneities. Thus, the photometric observations 
were acquired over a short time span to minimise possible 
variations in the out-of-eclipse shape caused by starspot 
migration. Also, special care was taken to schedule the 
observations near the eclipses so that a very dense phase 
coverage could be achieved. To compute orbital phases, 
we adopted the ephemeris of D99 but correcting the spec¬ 
troscopic reference epoch to the primary eclipse: 

T(Min I) = HJD2450208.5068 + 2.771468 E 

The phased light curves are presented in Fig. [[]. Note 
the dense phase coverage, the relative shallowness of the 
eclipses and the out-of-eclipse variations presumably due 
to starspots. The individual photometric observations in 
the R (n = 2157) and I (n = 2188) bands are provided 
in Tables 1 and 2 (HJD, phase and Arnag), respectively, 
which are only available in electronic form. 

As mentioned above, CU Cnc has an optical compan¬ 
ion at ~12". The diaphragm used for the photometric ob¬ 
servations was 45" in diameter, so the light of the com¬ 
panion was included in the measurements. To estimate 
exactly the fraction of contributed light we made use of 
CCD images taken with the ALFOSC instrument at the 
Nordic Optical Telescope (La Palma). Several images in 
the V and I passbands were kindly acquired by R. J. 
Irgens on May 10, 1999 (HJD 2451309.375) at an orbital 
phase of 0.215. The reduction was carried out in the clas¬ 
sical manner and aperture photometry was performed to 
estimate the differential magnitude between CU Cnc and 
its companion. We obtained values of AU = 1.45 ± 0.01 
and A I = 1.20 ± 0.01, which translates into the following 
fractions of 3rd light, that can be directly incorporated 
into the analysis of the light curves: 

F 3 v = 0.21 ± 0.01; F 3 J = 0.25 ± 0.01 

From these we estimate the fraction of 3rd light in the R 
passband to be F jf* = 0.23 ± 0.01. 


3. Light curve analysis 

The fits to the light curves were performed using an im¬ 
proved version of the Wilson-Devinney program (Wilson 
& Devinney 1971; hereafter W-D) that includes a model 
atmosphere routine developed by Milone et al. ( |1992| ) for 
the computation of the stellar radiative parameters. A de¬ 
tached configuration was chosen for all solutions. Both 
reflection and proximity effects were taken into account, 
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Fig. l. R and / light curves of CU Cnc. A zero point shift of —1.45 has been applied to the R photometry for display 
purposes. The solid line in the bottom panel is the best-fitting synthetic light curve generated by W-D. The top panels 
show the residuals of the fits. Enlargements of the fits to the eclipse phases are shown in Fig. |j. 


even though the light curves do not show strong evidence 
for these. The bolometric albedo was set at the canoni¬ 
cal value of 0.5 for stars with convective envelopes. The 
gravity-brightening coefficient was adopted as 0.2 follow¬ 
ing the theoretical results of Claret (2000a). A mass ratio 
of q = Mb/Ma = 0.9184 was adopted from D99 and 
the temperature of the primary component (eclipsed at 
phase 0.0) was set to 3160 K (see Sect. [Ijfor discussion). 
For the limb darkening we used a logar ithmic law as de¬ 
fined in Klinglesmith & Sobieski (1970). In our W-D im¬ 
plementation, first- and second-order limb darkening co¬ 
efficients are interpolated at each iteration from a set of 
tables computed in advance using a grid of Kurucz model 
atmospheres. However, CU Cnc has a temperature some¬ 
what below the coolest model and a extrapolation had 
to be made. We also ran some tests by fixing the coeffi¬ 


cients at the values computed by Claret (2000b) for the 
appropriate effective temperatures. The change had negli¬ 
gible effects but we finally adopted this latter prescription 
to avoid extrapolations. Finally, third light fractions were 
set to the values provided in Sect. 0 . 


The iterations with the W-D code were carried out 
automatically until convergence, and a solution was de¬ 
fined as the set of parameters for which the differential 
corrections suggested by the program were smaller than 
the internal probable errors on three consecutive itera¬ 
tions. As a general rule, several runs with different start¬ 
ing parameters are used to make realistic estimates of the 
uncertainties and to test the uniqueness of the solution. 


Simultaneous W-D fits to the RI light curves were 
carried out using the mean error of a single measure¬ 
ment as the relative weight, w\, of each passband (“curve- 
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dependent” weighting scheme). Note that the number of 
photometric measurements we have is so large that 4-point 
normals had to be computed at the densest phases in or 
around the eclipses when running the W-D fits. We ini¬ 
tially solved for the following light curve parameters: the 
orbital inclination (i), the temperature of the secondary 
(T e ff B ), the gravitational potentials (Da and fie), the lu¬ 
minosity of the primary at each passband (La), a phase 
offset (A<^>), and the spot parameters. The orbital eccen¬ 
tricity was set to zero as the light curves do not show 
any evidence for eclipses of different width or a secondary 
eclipse at an orbital phase other than 0.5. Some tests were 
run to check for the possibility of a small non-zero orbital 
eccentricity. The W-D solutions converged towards a value 
of e = 0.001 ± 0.003, thus consistent with our adoption of 
a null eccentricity. 

As expected for a system with partial eclipses and sim¬ 
ilar components, numerous test solutions revealed that the 
ratio of the radii (k = tbIta, where and re are the 
fractional radii in units of the separation) is poorly con¬ 
strained. The problem is, in principle, severe because it 
implies that the light curves are not sufficiently sensitive to 
discriminate between the sizes of the components^. Thus, 
some source of external information has to be used to infer 
a value for k. With the ratio of temperatures well deter¬ 
mined from the light curves, a spectroscopic estimation of 
the luminosity ratio would provide the necessary informa¬ 
tion. 

Two high resolution and S/N echelle spectra were ob¬ 
tained with the UES instrument at the William Herschel 
Telescope (La Palma) as part of the Service Programme. 
The observations were carried on March 8-9, 2001 (HJD 
2451977.499) at orbital phase 0.287. The raw images were 
reduced using standard NOAO/IRAF tasks (including 
bias subtraction, flat field correction, sky-background sub¬ 
traction, cosmic ray removal, extraction of the orders, 
dispersion correction, merging, and continuum normaliza¬ 
tion). Equivalent width ratios for the components were 
measured for 10 relatively clean and isolated Ca I and 
Fe I absorption features. The average of the individual 
values yielded EWb/EWa = 0.73 ± 0.05 at a mean wave¬ 
length of A = 5900 A. Since the effective temperatures of 
both components are very similar, as shown below, the 
equivalent width ratio is analogous to luminosity ratio. 
From this we derive values of [L b /La]v = 0.72 ± 0.05, 
[L b /L a ]r = 0.74 ± 0.05 and [L b /L a \i = 0.76 ± 0.05, for 
the luminosity ratios at the V, R and I passbands, respec¬ 
tively. Our value for the U-band luminosity ratio is indeed 
in good agreement with that obtained by D99. Although 
not explicitly given in the paper, Delfosse (2000, priv. 
comm.) kindly computed the luminosity ratio from their 
spectra and obtained a value of 0.68, which lies within one 
sigma of our determination. 


1 Note, however, that the sum of fractional radii is extremely 
well determined because it depends strongly on the total du¬ 
ration of the eclipses. 


In practice, the luminosity ratio or the ratio of radii 
is not one of the parameters considered explicitly in the 
W-D code. Instead, are the surface potentials of the stars 
(which, along with q, control the sizes) that can be fixed 
or left free. We thus ran a number of solutions by adopting 
different values for the potential of the primary star until 
reaching agreement between the resulting luminosity ratio 
and the observed one. This occurred for Da = 19.0. All 
further solutions were ran by fixing Da at this value. 

The shape of the out-of-eclipse region in the light curve 
clearly indicates that one or both components have sur¬ 
face inhomogeneities. A feature incorporated by the W-D 
program is the capability of fitting a simple spot model. 
The spots, assumed to be circular in shape, are described 
by W-D through four parameters, namely, the longitude, 
latitude, angular radius, and temperature ratio (spot rel¬ 
ative to photosphere). Some of the spot parameters must 
be fixed when running solutions to prevent lack of con¬ 
vergence. For example, the latitude of the spots can be 
inferred theoretically. Granzer et al. ( 200C| ) find that a 
star such as the components of CU Cnc (0.4-Mg star 
in the ZAMS with an angular velocity 10 times the so¬ 
lar value) favours strongly the formation of spots at lati¬ 
tudes of about 60°. Tests ran with W-D indeed indicated 
a marginally better fit for spots at high latitudes. 

As it seems obvious, the radius of the spot and the 
temperature contrast are strongly correlated. No empiri¬ 
cal information is available on the temperature ratio be¬ 
tween the spots and the surrounding photosphere for CU 
Cnc. However, Hatzes (1995) analysed Doppler tomogra¬ 
phy of the M-type eclipsing binary YY Gem and found the 
presence of dark areas cooler than the photosphere with 
Tspot/Tphot = 0.87 — 0.92. For CU Cnc, we expect simi¬ 
lar values. We thus ran a series of solutions and tested a 
variety of spot combinations (number of spots, locations, 
sizes, temperature ratios, etc) to achieve the best possible 
fit to the light curves (see a detailed discussion of a sim¬ 
ilar procedure in TR02). Our preferred solution calls for 
two spots - one on each component -. The spot on the 
primary component, which transits the central meridian 
at orbital phase 0.438, appears to be relatively small, with 
a radius of 9°, and ~450 K cooler than the photosphere. 
In contrast, the secondary component has a much larger 
spot with a radius of 31° and a temperature difference 
with the surrounding photosphere of ~200 K. The center 
of the dark area transits the central meridian at orbital 
phase 0.997. Its large size and relatively small tempera¬ 
ture difference seems to indicate that this is more likely 
a spot complex rather than a homogeneous spot. Thus, 
the parameters determined represent an average over an 
extended photospheric area probably covered with patchy 
dark spots. 

The photometric effect of the spot configuration 
adopted here is illustrated in Fig. ||. Note that the spot 
solution has only a mild effect on the rest of the light 
curve parameters. Furthermore, the spot solution is not 
unique and other configurations may be equally valid from 
a strictly numerical point of view. 
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Fig. 2. Photometric effects of the spots on the R and I 
light curves. The spike during the phase of the primary 
minimum arises when the primary component is partially 
eclipsed and the spot on the secondary component is in 
front thus enhancing its relative contribution. 

The best-fitting two-spot model described above yields 
rms residuals of 0.008 mag and 0.005 mag for the R and I 
light curves, respectively. These small values reflect both 
the quality of the photometry and the excellent perfor¬ 
mance of the spot model. The synthetic light curve super¬ 
imposed to the observations can be seen in Fig. |l], together 
with the residuals of the fit for each light curve. Also, for 
illustrative purposes, we show a detail of the eclipse phases 
in Fig. |j. The solution indicates that the eclipses are par¬ 
tial and that ~25% of each component’s flux is blocked 
during the corresponding eclipses. The parameters result¬ 
ing from the light curve fit are listed in Table |^. The un¬ 
certainties given in this table were not adopted from the 
formal probable errors provided by the W-D code, but 
instead from numerical simulations and other consider¬ 
ations. Several sets of starting parameters were tried in 
order to explore the full extent of the parameter space. 
In addition, the W-D iterations were not stopped after a 
solution was found, instead, the program was kept run¬ 
ning to test the stability of the solution and the geome¬ 
try of the x 2 function near the minimum. The scatter in 
the resulting parameters from numerous additional solu¬ 
tions yielded estimated uncertainties that we consider to 
be more realistic, and are generally several times larger 
than the internal statistical errors. For CU Cnc, however, 
an additional error source must be considered and this is 
the uncertainty in the measurement of the luminosity ra¬ 
tio. The empirical value obtained from the spectra carries 
a small but significant measurement error that must be 
properly accounted for. We did this by running new W-D 
solutions for values of set to yield the observed lumi¬ 
nosity ratio plus and minus 1 a. The average differences 
between the adopted and the new parameters were com¬ 
bined quadratically with the uncertainties described above 
to yield the values listed in Table |^. As a consequence of 
the uncertainty in the luminosity ratio, the radii of the 
components have larger errors (^2%) than would be ex¬ 
pected from the quality of the light curve 
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Fig. 3. Details of the fits to the primary and secondary 
eclipses. Note the high density of observations and the 
excellent agreement with the synthetic light curve. 

The only independent determination of the light curve 
parameters of CU Cnc comes from D99. As mentioned 
earlier, the authors only had sparse observations and the 
quality of the solution is not very high. Importantly, the 
effect of star spots were taken into account in the analysis. 
The final radii published by D99 have errors of 12-16%. 
Even so, the agreement between D99’s solution and the 
parameters listed in Table || is mostly within their large 
quoted uncertainties. The most significant difference is in 
the relative radii of the components, where D99 report 
very unequal components (fc = 0.67). Recall that the ratio 
of radii cannot be determined from the light curve alone 
so it is not surprising that D99 found rather inconsistent 
values. The sum of the relative radii, however, is tightly 
constrained by the duration of the eclipses and our value 
(0.1054±0.0006) is extremely close to that found by D99 
(0.1052). 

4. Absolute dimensions and effective 
temperatures 

The absolute dimensions of the CU Cnc components can 
be readily determined by combining the photometric pa¬ 
rameters in Table ^ with the spectroscopic solution from 
D99. When doing so, we obtain the physical properties 
listed in Table [| 

The absolute value of the effective temperature for 
eclipsing binaries cannot be determined from the light 
curves alone and independent methods must be used. 
Temperature determination is usually quite problematic 
when one has to deal with very cool stars such as CU 
Cnc. Two independent approaches present themselves to 
estimate the effective temperature of CU Cnc: one based 
on photometric calibrations and one based on the known 
absolute dimensions and distance. 
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Table 3. Results from the light curve analysis for CU 
Cnc. 


Parameter 

Value 

Geometric and radiative parameters 

P (days) (fixed) 

2.771468 

e (fixed) 

0 

i (deg) 

86.34A0.03 

q = Mb /Ma (fixed) 

0.9184 

12a (fixed) 

19.00A0.22 

12b 

19.39A0.40 

va (volume) 

0.0553±0.0007 

re (volume) 

0.0501A0.0012 

k = r B /r A 

0.906A0.025 

±eff (K) (fixed) 

3160 

ri-iB /ctuA 
eff / -*■ eff 

0.9892A0.0019 

L b /La ( R band, phase 0.287) 

0.74A0.05 

Lb/La ( I band, phase 0.287) 

0.76±0.05 

±3 ( R band, phase 0.215) (fixed) 

0.23 

±3 (/ band, phase 0.215) (fixed) 

0.25 

Albedo (fixed) 

0.5 

Gravity brightening (fixed) 

0.2 

Limb darkening coefficients (Logarithmic law) 

xa and i/a ( R band) 

0.875 , 0.390 

x B and y B (R band) 

0.875 , 0.390 

xa and i/a (I band) 

0.844 , 0.530 

xb and j/b ( I band) 

0.844,0.530 

Spot parameters 

Phase for spot # 1 

0.438 

Radius for spot #1 (deg) 

9 

AT for spot #1 (K) 

450 

Location of spot #1 

Primary 

Phase for spot # 2 

0.997 

Radius for spot #2 (deg) 

31 

AT for spot # 2 (K) 

200 

Location of spot # 2 

Secondary 

r.m.s. residuals from the fits 

cm (mag) 

0.008 

a: (mag) 

0.005 


Synthetic colors, which are the base of photometric 
calibrations, are computed by modelling the stellar atmo¬ 
spheres. Albeit some discrepancies remain, M-star model 
atmospheres have made significant progress over the past 


Table 4. Absolute dimensions and radiative properties for 
the components of CU Cnc. 


Property 

Primary 

Secondary 

Spectral type 

M3.5 Ve 

M3.5 Ve 

Mass (Mg) 

0.4333A0.0017 

0.3980A0.0014 

Radius (Rq) 

0.4317A0.0052 

0.3908A0.0094 

log fit (cgs) 

4.804±0.011 

4.854A0.021 

Teff (K) 

3160A150 

3125A150 

log(L/L 0 ) 

— 1.778A0.083 

-1.884A0.086 

A/boi (mag) 

9.19A0.21 

9.45A0.21 

7r (mas) 

78.05A5.69 

A/v (mag) 

ll.95A0.16 

12.31A0.16 

BCy (mag) 

-2.76A0.26 

—2.86A0.26 


decade (Allard et al. 1997 ). Line opacities from a large 
number of molecular species are now included in the cal¬ 
culations and substantial improvements have been made 
in the input physics, particularly the equation of state. 
For cool stars, the infrared colors are especially suited to 
estimate the effective temperature because most of the 
stellar flux is emitted in that wavelength range. Thus, we 
collected from the literature the available magnitudes and 
colors for CU Cnc, which are listed in Table |[ Simple 
linear transformations were applied to refer the observed 
magnitudes and colors to the photometric systems most 
widely used. 

The magnitudes and color indices in Table corre¬ 
spond to the combined light of the two components of CU 
Cnc. Since the effective temperatures of the two compo¬ 
nents are very similar, it is justified to use the joint color 
indices as those representing the “mean” component of 
CU Cnc. To estimate the mean effective temperature we 
employed a number of modern photometric calibrations. 
We assumed no interstellar reddening in our calculations 
as expected for a system that lies less than 20 pc away. 
The resulting temperatures, together with the references, 
atmosphere model names, and comments, are presented 
in Table [| Most of the calibrations are based upon syn¬ 
thetic colors and model atmosphere calculations, except 
for Leggett (1992), who provides an empirical calibration. 
There is remarkably good agreement between all the tem¬ 
perature estimates. A plain average of the independent 
values (i.e., excluding one of the redundant temperatures 
from the BaSeL model) yields T e g = 3140 ±40 K. The for¬ 
mal error, however, does not reflect any possible systemat- 
ics that could be present in the calibration. We have thus 
adopted a more realistic error value of T e g = 3140 ± 150 
K based on the results of Leggett et al. ( 1996| ). 

As mentioned above, there is an alternative approach 
to compute the temperature of the com ponen ts based on 
the method described in Ribas et al. ( |l998| ). This is a 
simple method that relates the observed radius and lu¬ 
minosity with the temperature. In the case of CU Cnc, 
the absolute radii of the components are given in Table 
y]. The empirical luminosities were computed from the ap¬ 
parent mag nitud es, the parallax in the Hipparcos cata¬ 
logue (ESA 1997 ) and a bolometric correction. The latter 
is especially delicate because bolometric correction cali¬ 
brations are known to be the source of systematic errors. 
To minimise these, we used the bolometric correction in 
the K band, which was computed from the models of 
Allard et al. ( [2000| ) and resulted in BCkbb = ±2.75±0.08 
mag. This leads to an apparent bolometric luminosity of 
m bo i = 9.42 ± 0.08. From the luminosity ratio, the indi¬ 
vidual bolometric magnitudes can be computed and we 
find TOboiA = 10.04 ± 0.08 and mboiB = 10.31 ± 0.08. 
The absolute bolometric magnitudes follow from intro¬ 
ducing the Hipparcos distance: MboiA = 9-50 ± 0.18 and 
A^boiB = 9.77 ± 0.18. The effective temperature of each 
component can be computed in a straightforward manner 
and we obtain T e gA = 2950±120 K and T e g B = 2910±120 
K. 
















I. Ribas: The 0.4-Mg Eclipsing Binary CU Cnc 


7 


Table 5. Joint observed magnitudes and color indices for CU Cnc. 


Mag. or color 

Value 

V 

11.89 ±0.01 

(R -1) c 

1.60 

( y-i)b 

2.80 

JW ASS 

7.536 ±0.018 

H 2MASS 

6.912 ±0.030 

Admass 

6.626 ± 0.029 

A'bb 

6.670 ± 0.031 

A cit 

6.650 ± 0.031 

(U - I<)b b 

5.22 ±0.03 

(V — A)c it 

5.24 ±0.03 


Notes 


From Weis (L991) 


From Weis (1991) and transformations in Bessell (1979) 
From Weis (1991) and transformations in Bessell (1979) 
2MASS Second Incremental Data Release b 
2MASS Second Incremental Data Release* 3 
2MASS Second Incremental Data Release 13 


Using transformations from Carpenter (2001) 
Using transformations from Carpenter (2001) 


Cousins (1976) photometric system 


http://www.ipac.caltech.edu/2mass/releases/second/doc/explsup.html 


c Bessell & Brett 
d Caltech photometric system (Elias et al. 1982, 1983) 


homogenized system 


Table 6. Effective temperature determinations for the mean component of CU Cnc. 


Reference 


Allard et al. (2000] 


20001 ) 

al. (|199£ 


Leggett (1992) 


Hauschildt et al. ([1999 ) 
Lejeune et al. (1998) 
Lejeune et al. (1998) 
Bessell et al.(1998) 


Model T eff (K) 

STARDusty2000 3130 

NextGen 3170 

BaSeL 3.1 3170 

BaSeL 2.2 3170 

NMARCS 3070 

3145 


Colors used 
(V - I)c, (V - JOcit 

(y — a')cit 

(V - I)c, (V - K )BE 
0 V - /)c, (V - K )BB 
(V - I)c, (V - K) BE 
(R-I)c, (V-I)c, (V-K) cit 


The agreement with the photometric determinations 
described above is not perfect but nevertheless fairly good 
given the uncertainties. The values agree within one sigma 
of their respective error bars. However, a caveat is due at 
this point. As it turns out, there are reasonable concerns 
about the accuracy of the Hipparcos parallax, even taking 
into account its relatively large error (see Table |]). CU 
Cnc is a faint object near the threshold for detection by 
Hipparcos. Thus, the individual transit data are of poor 
quality. A reanalysis of the transit data (considering also 
the visual component) has not provided any better results 
(Arenou 1999, priv. comm.). As a reference, we have com¬ 
puted the temperature of CU Cnc as a function of the 
distance^: 

T effA = 2610^^ K; T effB = 2575^^ K 

Given the concerns with the Hipparcos distance raised 
above, we prefer to adopt the temperatures derived from 
photometric calibrations. The final adopted individual 
temperatures for CU Cnc (computed from the average 
temperature and the temperature ratio derived from the 
light curves) are listed in Table ^ along with their asso¬ 
ciated errors. Also included in Table ^ are the absolute 
magnitudes of the components in the V band, which have 
been computed from the parallax listed in the Hipparcos 

2 Note that there is a weak dependence of the bolometric 
correction on the temperature, so these expressions are only 
valid for temperatures not too far off from those in Table 


catalogue. Note that an empirical determination of the 
bolometric correction for each component follows from the 
comparison of the U-band absolute magnitude with the 
bolometric magnitudes (resulting from the radii, effective 
temperatures and Mboi® = 4.74). 

5. Other observational information 

5.1. Age and chemical composition 

One of the goals of the present study is to carry out a 
critical analysis of low-mass stellar models similar to that 
in TR02. Thus, to make the test as stringent as possi¬ 
ble, one has to keep the number of degrees of freedom at 
a minimum. Two of the key parameters that determine 
the observable physical properties of a star are its age 
and chemical composition (to simplify, metallicity). The 
metallicity of a star is generally determined through high- 
resolution spectroscopy and comparison with synthetic 
spectra. Unfortunately, this is virtually impossible for M 
stars because of the huge quantity of spectral features and 
the shortcomings of the atmosphere models. In turn, the 
age (and eventually the metallicity) can only be derived 
if the star belongs in a well-studied cluster. Field M stars 
are, thus, not ideally suited for reliable estimations of age 
and metallicity. 

One possible way to obtain a rough estimation of the 
age and chemical composition of a field M star such as CU 
Cnc is through the analysis of its space velocities. High ve¬ 
locity would likely indicate an (old) halo population and 
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low velocity a (young) star in the galactic disk. The space 
velocities (U, V, W) of a star are readily computed from its 
position, proper motions, radial velocity and distance. In 
the case of CU Cnc, we have adopted the Hipparcos posi¬ 
tion, proper motion and distance, and the radial velocity 
has been taken from D99. The resulting space velocities 
are U = —10.7 km s _1 , V = —4.9 km s -1 , and W = —10.6 
km s -1 . Note that we have followed the convention where 
positive values of U , V, and W indicate velocities towards 
the galactic center, galactic rotation and galactic North 
Pole, respectively. 

The low values of these velocities are an indication 
of a disk population and thus a relatively young star. 
A closer look at the values revealed that these are strik¬ 
ingly similar to the v elociti es of the Castor movi ng gro up 
(see Anosova & Orlov 1991 Barrado y Navascues 1998|) o f 
which YY Gem is a likely member. Ribas & Rebolo ( [2002 ) 
utilised various criteria (isochrones, rotation-activity rela¬ 
tionships, etc) to estimate an age of 320±80 Myr and a 
metal content of Z = Z Q for the moving group mem¬ 
bers. The mean space velocities of the moving group have 
been recently re-determined by Ribas & Rebolo ([2002) to 


in clusters (e.g., Duncan 1981; Soderblom et al. 


km s -1 , and < W >= —9.4 ± 2.1 km s -1 . These are in¬ 
deed very similar to the values found for CU Cnc, which 
strongly suggests its membership of the Castor moving 
group. Thus, it is sensible to assume that CU Cnc has 
the same age and chemical composition as the rest of 
the group members - similarly to what is routinely done 
in stellar cluster studies. As it happened with YY Gem 
(TR02) the fortunate fact that CU Cnc belongs in the 
moving group and has a well-determined age and metal 
content greatly enhances the value of the star since it 
allows for an unusually stringent test of the theoretical 
models. 


5.2. Lithium abundance 

As it has been known for some time, the measured abun¬ 
dance of lithium in the stellar atmosphere can be used 
as an age indicator for young cool stars, parti cularl y 

1990|) . 


pears to have been mostly destroyed, yet not completely. 
This is in apparent contradiction with Li destruction se- 
quen ces in clusters and as sociat ions (Barrado y Navascues 
et al. 1999; Stauffer et al. 1999), which indicate that mid- 
M type stars fully deplete their initial Li abundance in as 
little as a few times 10' years. (Recall that the estimated 
age for CU Cnc is ~320 Myr.) Barrado y Navascues et al. 
(1997) observed a similar discrepancy with the observed Li 
abundance (log./V(Li) = 0.0) for the eclipsing binary YY 
Gem. In this case, Li also appears to be severely depleted 
but not completely destroyed, as would be expected for 
YY Gem’s age of ~320 Myr. If we rely on the ages deter¬ 
mined for these binaries, the two results seem to indicate 
that Li depletion for M stars in binary systems might be 
not as efficient as for single stars. For example, orbital 
synchronization due to tidal forces could be responsible 
for inhibiting turbulent mixing, thus slowing down Li de¬ 
pletion (Barrado y Navascues et al. 1997 ). The tentative 
detection of Li in the spectrum of CU Cnc is a puzzle yet 
to be resolved a detailed discussion of this topic is left for 
an upcoming work. 


5.3. Stellar activity 

Other observational information for CU Cnc can be in¬ 
ferred from the high-resolution spectrum described in 
Sect. ||. An inspection of the spectrum reveals very strong 
Ha and H/3 emission features (double lines). The measure¬ 
ments yield (continuum-corrected) Ha equivalent widths 
of 3.85 A and 4.05 A for the primary and secondary com¬ 
ponents, respectively. For comparison, Young et al. ( 1989 ) 
report an excess equivalent width of ss2 A for the presum¬ 
ably coeval but more massive system YY Gem. The Ha 
equivalent widths of both CU Cnc and YY Gem are rather 
large but not unreasonable when compared with young M- 


Theoretical studies demonstrate that the 'Li nuclei are 
destroyed by proto n colli sions at temperatures a bove ~ 2.5 
MK (Bodenheimer 1965 ; D’Antona & Mazzitelli 1984 ). In 
low mass stars, mixing through convection processes is 
so efficient that Li is depleted at a very rapid pace. We 
have inspected our high-resolution spectrum described in 
Sect. U and we report a tentative detection of weak Li I 
A6708 features at the expected wavelengths for both the 
primary and secondary components. The equivalent width 
measurement (corrected for light dilution) yields a value 
of ~50 mA for both stars. From this and the theoretical 


type stars (see, e.g., Soderblom et al. 1991). Note, how¬ 
ever, that the strong Ha emission in these binary systems 
is more related to the spin-up caused by orbital synchro¬ 
nization rather than age. 

In addition to displaying Ha chromospheric emission, 
CU Cnc is also both and EUV and X-ray source. From 
the observations in the ROSA T All-Sky Bright Source 
Catalogue (Voges et al. 199S) and the calibration in 


study of Pavlenko et al. (1995), one can roughly estimate 
the Li abundance to be log N (Li) ss —1.1 (in the scale 
in which log./V(H) = 12), and k>g[n(Li)/n 0 (Li)] = —4.1. 
Thus, the initial Li content of CU Cnc’s components ap¬ 


Schmitt et al. ( |1995| ) one obtains an integrated X-ray lu¬ 
minosity of logLx(erg s -1 ) = 29.03 for CU Cnc. Again, 
this value can be compared with YY Gem’s, which was 
estim ated to be log Lx(erg s _1 ) = 29.61 (Giidel et al. 
1993). One must keep in mind that the size of the stars 
in CU Cnc is significantly smaller than those of YY Gem. 
Thus, the ratio Lx/Aboi is a more realistic measure of 
the stellar activity or, in other words, of the efficiency 
in producing high energy emissions. The calculations of 
Lx/Lboi yield a value very close to 10~ 3 for both CU Cnc 
and YY Gem, thus indicating very similar activity levels. 
Numerous studies (see, e.g., Jeffries & Tolley 199S| ; Gagne 
et al. 1995) suggest that this value of Lx/Aboi = 10~ 3 
indeed marks a saturation limit in the stellar activity. 































I. Ribas: The 0.4-Mg Eclipsing Binary CU Cnc 


9 


6. Critical comparison with stellar models 

CU Cnc is only the third M-type eclipsing binary for which 
accurate absolute dimensions (masses and radii) have been 
determined. In addition, we have been able to set rather 
stringent constraints on the effective temperatures of the 
components as well as on the age and chemical compo¬ 
sition of the system. With this information in hand, we 
are in a position to carry out a critical evaluation of the 
stellar model predictions. For the first time we can pro¬ 
vide observational checks on M stars of about 0.4 Mq for 
which no reliable empirical information has been available 
to date. The procedure followed here is very similar to that 
described in TR02 for YY Gem, but with a significant dif¬ 
ference. While the two components of YY Gem are almost 
identical, this is not the case for CU Cnc where the compo¬ 
nents differ by almost 10% in mass. Since the two stars in 
the binary system are coeval, a single isochrone expected 
to fit both components simultaneously. This implies that 
not only the position in the observational diagrams but 
also the slopes of the isochrones come into play. Thus, the 
models must reproduce both the absolute and the relative 
location of the components. 

For our comparison we have considered nine differ¬ 
ent sets of theoretical c alcula tions: Swenso n et a l. ( 1994 ), 
D’Antona & Mazzitelli ( 1997 ), Siess et al. ( |l997 ), Baraffe 
et al. (199b), Palla & Stahler ( 1999 ), Charbonnel et 
al. ( 1999 ), Girardi et al. ( 2000|) , Yi et al. ( flOOl ), and 


some of the models yield reasonably good, albeit not per¬ 
fect, fits to the data. The three isochrones that exhibit the 


all models except those by Palla & Stahler (1999), with an 
oldest age available of 100 Myr (this has no effect whatso¬ 
ever on our conclusions). The models of Siess et al. (1997) 
and Baraffe et al. ( [1998 ) employ the most sophisticated 
physics (equation of state) and boundary conditions (non¬ 
grey model atmospheres), so, in principle, they would be 
expected to provide the best fit to the physical properties 
of CU Cnc. 

We have compared the observational data for CU Cnc 
with the model predictions at different planes. The plots 
are shown in the eight panels of Fig. ||. For clarit y rea - 
sons, we have plot ted th e models of Siess et al. ( 1997 ) 
and Baraffe et al. ( |1998| ) separated from the rest in the 
right panels. Note that we have also included the observa¬ 
tional data on YY Gem AB and V818 Tau B (both from 
TR02), which will help to understand the comparison in 
a wider context. The best determined parameters of our 
analysis are the masses and the radii, which are purely em¬ 
pirical, whereas the effective temperatures and absolute 
magnitudes hinge to some extent on external calibrations. 
Thus, the mass-radius diagram in Fig. [|a constitutes the 
most reliable check on the models, also because the error 
bars of the measurements are very small. As can be seen, 


Bergbusch & VandenBerg (2001). These include essen¬ 
tially all of the modern models for low mass stars, even 
though some of them do not quite reach down to the 
masses of CU Cnc. We have, nevertheless, incorporated 
them so that their performance at higher masses can be 
compared with the rest of the models. Isochrones for an 
age of 320 Myr and Z = 0.018 have been interpolated for 


best perfor mance are those by Sw enson et al. ( 1994 ), Palla 
& Stahler ( 1999 ), and Yi et al. ( 2001 ). H oweve r, among 
these, only the models of Swenson et al. ( 1994 ) succeed 
in reproducing at the same time the relative location (i.e. 
slope) of the binary components. The models of Siess et 
al. (1997) and Baraffe et al. (1998), in spite of using state- 
of-the-art physical ingredients, do not yield a good fit to 
the data, although the latter perform slightly better. The 
radii that these models predict at the CU Cnc’s masses 
are about 10-14% smaller than observed. (Note that this 
difference amounts to ~7<r.) The plots indicate that the 
huge discrepancies between model predictions and obser¬ 
vations encountered by TR02 at slightly higher masses 
(YY Gem AB and V818 Tau B) are not quite as dramatic 
at the mass regime of CU Cnc, and some of the models 
produce reasonable fits to the data. 

The panels in Fig. [|b introduce the effective tem¬ 
perature in the comparison. Here the internal discrepan¬ 
cies among the seven models grow larger as we move to¬ 
wards lower masses. At around 0.4 Mq differences of up 
to 20% exist in the effective temperatures predicted by 
the different isochrones. Interestingly, most of the models 
seem to indicate temperatures that are significantly (10- 
15%) higher than observed for CU Cnc. The isochrones of 
Swenson et al. (1994) and Yi et al. (2001) clearly stand 
out among all and appear to produce a fairly close fit. But 
again, only the isochrone computed from the models by 
Swenson et al. (1994) reproduces the relative location of 
the components to achieve an outstanding fit to CU Cnc 
and also YY Gem, although V818 Tau B seems signifi¬ 
cantly cooler than model predictions. 

The mass-luminosity diagram in Fig. ^c: presents a sit¬ 
uation which is analogous to that described above. Most 
of the isochrones predict Afy’s for CU Cnc that are too 
bright by u p to 2 magnitudes. For ex ample , the models of 
Siess et al. ( 1997 ) and Baraffe et al. ( 1998 ) yield close fits 
to both V818 Tau B and YY Gem while huge differences 
are found with the observe d mag nitudes of CU Cnc. The 
models by Swenson et al. ( |1994| ) are able to both repro¬ 
duce the absolute magnitude of CU Cnc and YY Gem, 
while also achieving a quite close fit to V818 Tau B. At 
the same time, the slope of the isochrone matches perfectly 
the relative configuration of the CU Cnc components. 

Finally, the H-R diagrams in Fig. |]d do not reserve 
any surprises. Most of the models tend to overestimate 
the effective temperature at a given absolute magnitude. 
The best agreement is found for the models of Swenson 


et al. (1994) that produce an isochrone with the correct 
slope that reproduces well the observed parameters of CU 
Cnc and YY Gem. Note that the relatively close fit by 


the models of Baraffe et al. (1998) is fictitious because it 
occurs for masses that are far off from those determined 
for CU Cnc. 
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Fig. 4. The location of CU One’s components on several observational planes is compared with isochrones from nine 
different theoretical models (see legends and references in text). The age for the isochrones is 320 Myr and the 
metallicity adopted is Z = 0.018, except for the isochrones by Palla & Stahler ( 199ijj ) that are for an age of 100 Myr 
and solar abundance. In addition to the components of CU Cnc (represented with filled circles), we show the locations 
of YY Gem’s mean component at a mass of 0.6 Mq and the secondary component of V818 Tau with a mass of 0.76 Mg 
(both represented as open circles). 
























I. Ribas: The 0.4-Mg Eclipsing Binary CU Cnc 


11 


7. Is CU Cnc akin to other stars of similar mass? 


The comparison discussed above leaves no doubt that 
most of current stellar models fail in reproducing the ob¬ 
served physical pr operti es of CU Cnc. Only the models 
by Swenson et al. ( 1994 ) yield reasonably good (even ex¬ 
cellent) isochrone fits to the data. Even though CU Cnc 
is the only object with a mass around 0.4 Mg with fully 
characterised physical properties, there are other stars in 
this regime whose masses have been spectroscopically de¬ 
termined to good accuracy. This is the case of G1 570C, 
G1 623A, G1 644A, G1 644Ba, G1 661A, and G1 661B, all 
with masses between 0.34 and 0.42 Mg (see Forveille et al. 
199£; Segransan et al. 2000; Delfosse et al. 2000| ; Martin 


et al. 1998: Mazeh et al. 2001). We have compiled V and 


K absolute magnitudes for these stars from Delfosse et al. 
( 2000}) . Two mass-absolute magnitude plots are shown in 
Fig. |5|, one for the V magnitude and one for the K magni¬ 
tude. The results are quite striking as CU Cnc appears to 
be fainter in both plots than other stars of the same mass. 
The magnitude difference is about 1.4 mag in the V band 
and 0.35 mag in the K band. The apparent faintness of 
CU Cnc in the V b and ha d already been poi nted out by 
D99, Delfosse et al. ( 2000 ) and Mazeh et al. ( 2001 ). Note 
that t he isoc hrones computed with the models of Baraffe 
et al. ( 1998 ) and Siess et al. ( 1997 ) plotted in the figure 
produce a reasonably close fit to the stars other than CU 
Cnc in both observational diagrams. On the other hand, 
the isochrone by Swenson et al. ( 1994 ) - only available in 
My -, which closely matches CU Cnc, fails in reproducing 
the observed absolute magnitudes of the other stars. 

Let us carefully analyse possible reasons for the dis¬ 
crepancy. One scenario to explain the disagreement is the 
effect of enhanced stellar activity in CU Cnc due to bi- 
narity (higher rotational velocity). This would cause the 
appearance of large spots on the surface of the stars and 
also an increase of the TiO absorption. It is well known 
that the V passband is strongly affected by TiO absorp¬ 
tion whereas the K band is mostly free from TiO bands. 
This might cause the star to appear especially faint in 
the V band. There are, however, evidences that argue 
against this scenario. All observational data (X-ray lu¬ 
minosity, Ha emission) indicate that CU Cnc is not espe¬ 
cially active when compared with YY Gem, which is quite 
closely matched by the models. At CU Cnc’s mass range, 
the ROS AT All-Sky Bright Source Catalogue (Voges et 
al. 1999) yields an estimated X-ray luminosity for the 
stars in the Gliese 644 system (plotted in Fig. ||) of 
logLx(erg s -1 ) = 29.14, thus indicating an activity level 
very similar to CU Cnc’s. In addition, the available multi¬ 
band photometry (see Table ^j) refutes the idea that the 
apparent faintness of CU Cnc is only restricted to the V 
band. The spectral energy distribution deduced from the 
photometric data is internally consistent, i.e. all photo¬ 
metric indices yield the same effective temperature. Thus, 
it is the temperature of CU Cnc and not the photome¬ 
try the source of the discrepancy. Using (V — K) indices 
and the calibrations in Sect. 0 , we estimate that CU Cnc 



0.6 0.5 0.4 0.3 

Mass (M o ) 


Fig. 5. Mass-absolute magnitude diagrams in the V 
and K bands for the lower main sequence. CU Cnc 
is shown together with YY Gem and other stars with 
spectroscopically-determined masses (see text). The lines 
represent isochrones for the estimated age of CU Cnc com¬ 
puted with different models. 


is some 300 K cooler than other stars of about 0.4 Mg. 
That the mean photospheric temperature is reduced by as 
much as 10% due to enhanced stellar activity in CU Cnc 
seems, in principle, quite unlikely. More so if we take into 
account that stellar activity produces photospheric plages 
in addition to dark spots that would (partly or totally) 
compensate for the temperature reduction. 

The distance to CU Cnc is another possible source 
for discrepancies since, as we have discussed in Sect. [| it 
might be subject to systematic errors. However, a magni¬ 
tude differential in the distance modulus of the star would 
be constant over wavelength and affect identically both V 
and K bands. This is a strong argument against this sce¬ 
nario because we observe a much larger deviation in My 
than in Mr. Nonetheless, a problem with the distance 
could still be responsible for a fraction of the discrepancy. 

Our estimate for the metal abundance for CU Cnc 
comes from its likely membership of the Castor moving 
group. The possibility for a different chemical composi¬ 
tion can be not completely discarded because the star 
might not belong in the group after all or have different 
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abundances than other group members. A very high metal 
abundance of [Fe/H] = +0.5 was put forward by D99 as 
a possible explanation for the anomalous location of CU 
Cnc in the mass-absolute magnitude diagram. This was 
proposed on the basis of the isochrone behaviour in the 
mass-My diagram for metal abundances of [Fe/H] = 0 
and [Fe/H] = —0.5. We have carried out this analysis 
again with the new parameters for CU Cnc and the mod¬ 
els of Baraffe et al. (1998). If, as assumed by D99, the 
absolute magnitude changes are supposed to scale linearly 
with metallicity, a value even higher than +0.5 is necessary 
since the V magnitude difference between [Fe/H] = 0 and 
[Fe/H] = —0.5 at a mass of 0.4 Mq is only 0.8 mag. We 
recall that CU Cnc appears to be dimmer by as much as 
1.4 mag. Furthermore, this scenario does not explain the 
discrepancy in the K band because the different metal¬ 
licity isochrones run virtually on top of one another. An 
additional argument against the high-metallicity hypoth¬ 
esis is the lack of other stars in the solar neighbourhood 
with such high metal abundance. In the thorough study 
of Edvardsson et al. ( 1993 ) the highest metallicity found 
was [Fe/H ] = +0.24. The catalog of Cayrel de Strobel et 
al. ( [2001 ) lists only a handful of stars with [Fe/H] > +0.4 
and most of them appear to be misclassified. The high 
metallicity proposed to resolve the issue with the location 
of CU Cnc in the mass-My diagram seems, at this point, 
quite unlikely. 

A younger age than estimated for CU Cnc from its 
membership of the Castor moving group would neither 
explain the discrepancy. A low mass star becomes dimmer 
as it evolves from the pre-nrain sequence phase towards the 
ZAMS, where the minimum luminosity is reached. CU Cnc 
is already fainter than the other stars in the same mass 
range. Alternatively, it could be that the stars used for 
comparison with CU Cnc (G1 570, G1 623, G1 644, G1 661) 
were all much younger (~10-20 Myr) and contracting to¬ 
wards the ZAMS. Once again, this hypothetical situation 
is not in agreement with the observed mass-magnitude di¬ 
agrams because the magnitude change in the V and K 
bands are predicted by the models to be very similar (ac¬ 
tually, a somewhat larger magnitude variation in the K 
band). 

Thus far, there are no plausible scenarios to account 
for the observed magnitude differences between CU Cnc 
and similar mass stars. Even though it might appear as 
improbable, the presence of a certain amount of dust ab¬ 
sorption could reproduce the o bserved behaviour. Indeed, 
according to Fitzpatrick ( 1999 ), the relationship between 
Ay and Ak for interstellar dust is Ay ss 8 Ak- Thus, 1.4 
magnitudes of extinction in V would roughly correspond 
to 0.2 magnitudes in K. Such magnitude corrections would 
bring the two mass-magnitude plots into fair agreement. 
At a distance of only ~13 pc, it can be definitely ruled out 
that line-of-sight interstellar dust could account for over 1 
magnitude of extinction in the V band ( E(B — V) = 0.42). 
The alternative is therefore circumstellar dust. Disks of 
circumstellar dust have been directly detected in a num¬ 
ber of nearby stars (see, e.g., Greaves et al. 1998 and ref¬ 


erences therein) and some of them have estimated ages 
similar to CU Cnc’s. If the disk was coplanar with the 
binary’s orbit, as it might likely be assumed, we could be 
observing the stars through a line of sight with enhanced 
dust absorption. There are no observational grounds to 
support such hypothesis yet but, at this point, this is the 
only scenario that explains the observations. 

If the disk hypothesis is proven to be correct, our tem¬ 
perature estimate (which is based on multi-colour pho¬ 
tometry) and thus the comparison with stellar models in 
the temperature and absolute magnitude diagrams would 
have to be revised accordingly. A rough estimate made 
by correcting for the differences of 1.4 mag in V and 
0.35 mag in K raises the temperature of the components 
to T eSA = 3440 ± 150 K and T eSB = 3400 ± 150 K. 
Also, the effective temperatures obtained from the radii, 
the corrected apparent magnitudes and the Hipparcos 
parallax (see Sect. [|) are T e g A = 3290 ± 140 K and 
T e SB = 3240 + 140 K, thus achieving a remarkable 1 sigma 
agreement with the photometric values. 

8. Conclusions 

The lack of observational constraints is one of the most 
severe deficiencies that stellar models for the lower main 
sequence have to face. Eclipsing binaries are thus crucial 
because they yield simultaneous determinations of masses 
and radii (and also temperatures) upon which the stel¬ 
lar models can be checked. We have therefore acquired 
and analysed high-quality R and I light curves with very 
dense phase coverage (over 2100 observations) of the low- 
mass eclipsing binary CU Cnc. In the light curve analysis, 
carried out with the W-D program, we resolved an inher¬ 
ent indeterminacy of the ratio of radii of the components 
(caused by the partial eclipses) by fixing the luminosity 
ratio to the value obtained from a high-resolution spec¬ 
trum. Also, special care has been put in modelling the 
conspicuous out-of-eclipse variations through surface in- 
homogeities. The orbital and physical properties resulting 
from the light curves were combined with the spectro¬ 
scopic parameters from D99 to achieve precisions better 
than 2% in the absolute dimensions of CU Cnc compo¬ 
nents. As a result of this, CU Cnc is the third system 
(together with YY Gem and CM Dra) to reach sufficient 
accuracy in its absolute dimensions so that critical tests 
of stellar models can be performed. 

To make the comparison with models more stringent, 
we made use of further observational constraints. For ex¬ 
ample, the mean effective temperature of CU Cnc was 
determined from multi-band photometry (optical through 
infrared) and synthetic colours resulting from state-of-the- 
art atmosphere models. Additionally, the space motions of 
CU Cnc strongly suggest membership of the Castor mov¬ 
ing group. This allowed us to infer rough values for the 
age of the star (~320 Myr) and its chemical composition 
(Z ~ 0.018). With this host of observational information 
- quite uncommon for most eclipsing binaries -, we car¬ 
ried out a detailed comparison with all of the recent low- 
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mass stellar models. First, the most reliable check is in 
the mass-radius diagram because these are two quantities 
empirically determined from our analysis. The compari¬ 
son indicates that, whereas some of the models provide 
radius values that are in reasonable agreement with the 
observations, the most sophisticated models by Siess et al. 
(1997) and Baraffe et al. (1998) yield radii that are some 
10% smaller than observed. The study of three additional 
observational planes indicates that CU Cnc’s components 
are also significantly cooler (~15%) and fainter in the V 
band (~1.5 mag) than predicted by most models. 

A matter of some concern arises when comparing CU 
Cnc with other stars of similar mass. The comparison in¬ 
dicates that CU Cnc is some 1.4 mag fainter in the V band 
and 0.35 mag fainter in the I\ band. We have analysed a 
number of possible scenarios to explain the discrepancy 
and ruled out most of them on the basis of the available 
observational information. For example, we have discarded 
the effect of stellar activity, an error in the distance, a dif¬ 
ferent age or an extreme chemical composition. Only two 
scenarios appear to succeed in explaining the observed dif¬ 
ferences in the mass-magnitude diagrams: 1) CU Cnc’s 
components are, for a yet unknown reason, about 10% 
cooler than other stars of the same mass, and 2) there is a 
certain amount of (circumstellar) dust absorption that is 
responsible for the reddening of the spectral energy distri¬ 
bution. If the latter turns out to be the correct scenario, 
this is a very exciting result. Given the almost edge-on 
configuration of the eclipsing binary, we could be observ¬ 
ing CU Cnc through a coplanar dusty disk. Definite proof 
of this should come from submillimeter observations where 
the thermal emission of the dust could be detected. Note 
that CU Cnc belongs in the same moving group and prob¬ 
ably has si milar age to Vega and Fomalhaut (Barrado y 
Navascues 1998 ), both sta rs wit h bona-fide circumstellar 
dust disks (Holland et al. |l99<j ). With the tentative de¬ 
tection of a Vega-like edge-on dust disk and its relatively 
low mass and surface brightness, CU Cnc also becomes a 
very interesting target for planet searches through radial 
velocity studies or the transit method. 

In our comparative study, we have concluded that the 
evolutionary models of Swenson et al. ( 1994 ) clearly stand 
out among all because they achieve the best performance 
in reproducing all of the observed properties of CU Cnc 
and even the relative position of the components (much 
better constrained from the light curves). However, these 
models cannot obviously fit the position of other 0.4-Mq 
stars in the mass-absolute magnitude diagrams, whereas 
models such as Siess et al. ( 1997 ) and especially Baraffe 
et al. (1998) provide more reasonable fits. Unfortunately, 
until the discrepancy between the radiative properties of 
CU Cnc and stars of similar mass is resolved, no strong 
constraints on the models can be set in the temperature 
and absolute magnitude diagrams. A better determination 
of the parallax and an empirical estimation of the metal 
abundance would surely help to clarify the situation. 

Regardless of the issue with the temperature, the com¬ 
parison of CU Cnc’s absolute dimensions (empirically- 


measured masses and radii) with stellar model predictions 
clearly confirms the trend seen by TR02 for YY Gem and 
V818 Tau. The systematic underestimation of the radii of 
the stars can lead to severely underestimated stellar ages 
when using the H-R diagram diagnostics. In spite of re¬ 
cent advances, it seems clear that further adjustments are 
needed in evolution models for the lower main sequence to 
achieve a good description of the observed stellar physical 
properties. 
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